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We develop the general reconstruction scheme in two scalar model. The quintom-like theory which 
may describe (different) non-singular Little Rip or de Sitter cosmology is reconstructed, (fn)stability 
of such dark energy cosmologies as well as the flow to fixed points is studied. The stability of Little 
Rip universe which leads to dissolution of bound objects sometime in future indicates that no 
classical transition to de Sitter space occurs. 
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> 

O ■ I. INTRODUCTION 

^" 

' The observational data clearly indicate that current universe experiences the cosmic acceleration (dark energy 
. epoch). We are still far avifay from complete understanding of dark energy vifhich is often associated with some fluid 
with equation of state parameter w being close to —1, or with modified gravity. The observational data favor the 
ACDM model whose equation of state parameter w is equal to —1 exactly. Nevertheless, it is quite possible that the 
universe evolution is governed by phantom {w < —1) or quintessence (—1/3 > u; > —1) dark energy. If so then the 
' future universe can evolve to finite-time singularity like Big Rip [1] or quintessence-related soft singularity of one of 
three types according to classification [2] . The natural prescription to cure finite-time future singularity may be found 
(-H ' in frames of number of viable models of modified gravity [3]. However, it is still very interesting to understand if the 
I— 1 1 solution of singularity problem may be found within fluid dark energy. 

Recently, a new scenario to avoid future singularity has been proposed in [4] (for further development see Refs. [5- 
^ , 9]). According to this scenario the universe equation of state parameter w is less than —1, so that dark energy density 
■ increases with time but w sufficiently rapidly approaches —1 asymptotically. In this way, the finite-time singularity is 
QQ ' avoided. However, such proposed non-singular cosmology leads to dissolution of bound objects some when in future, 
Cf^ similarly to Big Rip singularity. That is why the scenario was called Little Rip cosmology. The scalar models to 
l/^ ' describe Little Rip were introduced in Ref. [5], they turn out to be phantom-like scalars. It is known that one scalar 
I , models are stable in phantom phase and instable in non-phantom phase. Moreover, the large instability occurs when 
T-H ' crossing the phantom-divide (or cosmological constant border). In order to understand Little Rip cosmology better, 
T-H , as well as its relation with asymptotically de Sitter universe and possible transitions between these two spaces, the 

' more realistic description of Little Rip may be necessary. 
^ I In the present paper we develop such description in terms of two-scalar tensor theory which represents kind of 
quintom model [10, 11] (for review, see [12]). Indeed, one of two scalars is taken to be phantom. 

In Section II, we consider a general formulation of reconstruction in two scalar model and investigate the stability 
of the solution. In this formulation, we construct a model which has a stable cosmological solution describing the 
phantom-divide crossing. In Section HI, we reconstruct a model which describes the cosmological solutions with 
and without Little Rip and investigate the (in)stability of the solutions. The existence of the solution describing de 
Sitter space-time and the stability of the de Sitter solution when it exists as well as possible transition of Little Rip 
cosmology to de Sitter one are investigated. In Section IV, we also consider the reconstruction of the two scalar model 
in terms of the e-foldings N and investigate the flow of the solution in terms of dimensionless variables, which give 
the fixed points for some solutions. Some summary and outlook are given in Discussion section. 
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II. RECONSTRUCTION OF TWO SCALAR MODEL AND (IN)STABILITY 

For the model with one scalar, the solution is stable in the phantom phase but unstable in the non-phantom phase. 
The instability becomes very large when crossing the cosmological constant line w = —1. In order to avoid this 
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problem, one may consider two scalar model [13, 14]. 
We now consider the following two scalar model 



S = J d^x^ { - l^i<l>)d^<l>d^ct> - \v{x)di.xd^X - V{<l>, x) | • (1) 
Here w(?!>) and r]{x) are functions of the scalar field (p and x, respectively. The FRW equations give 

uj{cP)^^+7j{x)x^ = -^H , V{cl>,x) = ^hH' + H) . (2) 



Then if 



ujit) + v{t) = - \nt) , Vit, t) = ^ {3f{tf + fit)) , (3) 



the explicit solution follows 

= X = ^ H = f{t). (4) 
One may choose that w should be always positive and rj be always negative, for example 

2 



vix) = -^^MxFTfW<o- (5) 

Here a is an arbitrary real function. We now define a new function x) by 

mx) = -y dM<l>) + 1 dxvix)^ , (6) 

which gives 

f{t,t) = m. (7) 

If ^(</',x) is given by using /((/), x) as 

the FRW and the following scalar field equations are also satisfied: 

= w{<P)4> + luj'{cl>)<p^ + SHcvicP)^ + ^^^^ , 

= V{X)X + lv'{x)x' + 3gr?(x)x + • (9) 

In case of the one scalar model, the instability becomes infinite at the crossing w = —1 point. One may expect that 
such a divergence of the instability does not occur for model with two scalars. 
By introducing the new quantities, X^, X^, and Y as 



X^^<A, X^^x, Y^l^, (10) 



the FRW equations and the scalar field equations (9) are rewritten as 

(X| - 1) - 3{X^ - Y) , 
{Xl-1)-3{X^-Y), 



dX^ 




dN 


2Huj{(l)) 


dX^ 


v'ix) 


dN 




dY _ 


3X^X^ ^1 


dN ~ 


Xdi + - 



H X^X^ + l^YjX^ + X^) 
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Here d/dN = H ^d/dt. In the solution (4), = X-^ = Y=1. The following perturbation may be considered 

X^ = 1 + SX^, X^ = 1 + 6X^, Y = 1 + 6Y . (12) 



Hence 




-^- 3 3 \ 

mI^x^ I , I -:^-3 3 

0-3-4/ 




(13) 



The eigenvalues of the matrix M are given by 



"'-w^-'- (») 



The eigenvalues (15) for the two scalar model are clearly finite. Hence, the instability, if any, could be finite and by 
choosing a in (5) properly, the instability can be removed, in general. In fact, right on the transition point where 
H = f'{t) = and therefore /'(^) = f'{x) = 0) for the choice in (5) with constant a, a{(j)) = a{x) = a > 0, we find 

i^{<t>) = -v{x)=^4 , ^'{<l>) = -^, l'ix)=0- (15) 



Then the eigenvalues (14) reduce to 



= ^ - 3 , = My = -3 . (16) 



Then as long as ^ < 3, all the eigenvalues are negative and therefore the solution (4) is stable. 

Hence, we gave general formulation of reconstruction in two scalar model and investigated the stability of the 
solution. By using this formulation, we can construct a model which has a stable cosmological solution corresponding 
to the phantom-divide crossing. 



III. RECONSTRUCTION OF LITTLE RIP COSMOLOGY 



In this section, by using the formulation of the previous section, we construct a model which may generate a Little 

Rip cosmology [4]. 

As an example we consider the following Hubble rate 

H = Hoe^' , (17) 

which corresponds to the Little Rip. Here Hq and A are positive constants. Eq. (17) shows that there is no curvature 
singularity for finite t. 

When we ignore the contribution from matter, the equation of state (EoS) parameter w of the dark energy can be 
expressed in terms of the Hubble rate H as 

2 FT 

- = -1-^. (18) 

Then if H > 0, w < -1. By using Eq. (18), one finds 

"^ = -1-^6-^*, (19) 

and therefore w < —1 and w — > — 1 when t — > +oo, and w is always less than —1 when H is positive. The parameter 
A in [4] corresponds to 2A/a/3 in (17) and is bounded as 2.74 x 10"^ Gyr"^ < A< 9.67 x 10"^ Gyr"^ by the results 
of the Supernova Cosmology Project [16]. 

In the model (17), H is always finite but increases exponentially, what generates the strong inertia! force. The 
inertial force becomes larger and larger and any bound object is ripped. This phenomenon is called a "Little Rip" [4] . 
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By choosing a in (5) as, 

a{t) = aoe-^* , (20) 
with a constant ao, we find and ri{x) in (5) as follows 

^i^) = ^ (V^"§ + ^'^0 - AHo) e^* , r,(x) = -^^«2+A2i/2e^x. (21) 

Using (6), we obtain 

m X) = j[- ^al + Xmi - XHo^ e^f + ^Ja^ + X^l e'^^ , (22) 
and the potential in (8) is given by 

n<P^x) =^\^^{- (V«o + A^//o'-AiJo) e^^ + ^Jal + X^Hie'^^' 

- [^Jal + A2£r2 - XH^ e^* + ^al + X'H^ e^^ | , (23) 

As written after (19), the parameter A is bounded as 2.74 x 10~^ Gyr~^ < 2A/\/3 < 9.67 x 10~^ Gyr~^ by the results 
of the Supernova Cosmology Project [16]. If we choose time so that the present universe corresponds to f = 0, we 
have Hq ~ 70km/s Mpc. 

For the model (17) with (20), the eigenvalues (14) are given by corresponding expressions when Xt 1, 



M^ = M^ = My = -i- Ae-^* , (24) 
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which arc negative. Therefore the solution is stable. 

Let us consider the possibility that the universe could evolve to the de Sitter space-time. In order for the solution 
corresponding to the dc Sitter space-time to exist, there should be an extremum in the potential and the potential 
should be positive there. If the extremum is local minimum with respect to (j), which is the canonical scalar and local 
maximum with respect to x, which is non-canonical or phantom scalar, the solution is stable. For the potential (23), 
there is an extremum when 



A' 
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■ + A2i/2 _ Ai/oj + ^al + X^H^ e^'^ = -y , (25) 
where the value of the potential ¥{(!), x) is given by 

mx) = -^, (26) 

which is negative and therefore there does not exist the solution corresponding to the de Sitter space-time. Hence, 
the universe does not evolve into the de Sitter space-time. 

We now show that the (asymptotically) Little Rip solution (17) is always asymptotically stable. For large t, one 
assumes the solution behaves as (17). In three eigenvalues (14), for the asymptotically Little Rip solution which is in 
the phantom phase H > 0, the eigenvalue My is negative. If we write a(x) in (5) as 

a{x) = XHoq{x)e^' , (27) 
the eigenvalue (14) can be expressed as a function of t = x as follows. 

If could be positive, q(t)q'{t) must be negative. Since q{t)'^ is positive, qit)"^ goes to a constant q{t)'^ ^ Qo > 0. 
Then due to the factor e'^* in the denominator of the first term in (28), the first term goes to small value for large t 
and we find ^ — 3 < 0. Therefore is asymptotically negative. For the eigenvalue M^, one gets 



V9(0"+i(V9(*)'+i-i) 



= ^-3- (29) 



5 



Then again in order that could be positive, q{t)q'{t) must be negative and therefore q{tY goes to a constant 
q{tY — Qo > 0- Due to the factor e'^* in the denominator of the first term in (29), the first term goes to small value 
for large t and — )• — 3 < 0. Therefore, all the eigenvalues are negative and the Little Rip solution is asymptotically 
stable. 

As another example, one can consider the following the model: 

H = Ha- Hie"^* . (30) 

Here Hq, Hi, and A are positive constants and we assume Hq > Hi and t > 0. Since the second term decreases when 
t increases, the universe goes to asymptotically de Sitter space-time. Then from Eq. (18), we find 

2Agie-^* 

w = -1 5-. (31) 

3{Ho-Hie-^tf 

As in the previous example (17), w < —1 and w ^ —1 when t +oo. In this model, there does not occur the Little 
Rip. The inertial force generating the Little Rip is given by 

Finer = mZo/a = ml(^H + H'^^ . (32) 

Here we consider two points separated by a distance I and assume there is a particle with mass m at each of the 
points. Since the magnitudes of H and H are bounded in the model (30), the Little Rip does not occur although the 

magnitudes of H and H become larger and larger in the model (17). 
For t 00, Eq. (31) gives the asymptotic behavior of w to be 



2XHie 



-xt 



---1-^^, (33) 

which is identical with (19) if we replace XHi/Hq with A. 
By choosing a in (5) as, 

a{t) = aoe"^* , (34) 
with a constant we find uji^)) and r](x) in (5) as follows 

w{cl>) = ^(j^7^!-XHi)e-^^, n{x) = -^Jal + Xm^e-^^. (35) 

Using (6), one gets 

Mx) = ifo - i | - (^^al + X^Hf-XHi^ e-^^ + + A^Ff e'^^^j , (36) 
and the potential (8) is given by 

mx)=4l4(^oA- 



al + X^Hl - XH^ e-^'t' - ^al + X^H^e'^^^ 
- (^^al + X^Hl - XH^ e-^"^ + ^ag + A2//2 g-^xj , (37) 



For the model (30) with (34), the eigenvalues (14) are given by 



M^ = M^ = -3+A, M^ = -3-^^i^. (38) 



Therefore as long as 3 > the solution (30) is stable. 

Since the solution can be unstable if 3 < , we again consider the possibility that the universe could evolve to the 
de Sitter space-time. For the potential (37), there is an extremum when 



(^^al + X^Hf - XHi^ e^^ + ^ al + X'^Hle>^^ = ~ Y + ^ 



(39) 
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where the value of the potential V{(l), x) is given by 

mx)= (-^ + Ai?o) . (40) 

If 12 > which is consistent with the condition 3 < that the solution (30) is stable, V{(f),x) is positive and 
there is a solution corresponding to the de Sitter space-time. Therefore there is a possibility that the universe could 
evolve into the de Sitter space-time. Note, however, the solution (39) corresponds to the minimum with respect to 
both of (j) and x and therefore the de Sitter solution is not stable. 
We may also consider the following the model: 

H = Hq + iJie"^* . (41) 

Here Hq, Hi, and A arc positive constants and we assume Hq > Hi and t > 0. Since the second term decreases when 
t increases, the universe goes to asymptotically de Sitter space-time. Then from Eq. (18), it follows 

2Affie-^* 

u; = -l+ 2- 42) 

3{Ho-Hie-^^f 

Hence, the EoS parameter is always larger than —1 and w ^ —1 when t — >■ +oo. Therefore the universe is in 
non-phantom phase. By choosing a in (5) as in (34), we find w(0) and r]{x) (5) as follows 

^('^) = ^ ( + >^'Hf + XHi^ e-^* , rj{x) = -^^al + Xm^e-^^. (43) 

Using (6), one gets 

f{<P,x) =Ho-jl^- {jal + \^Hl + \H^ e-^^ + + A^Ff g-Axj , (44) 
and the potential (8) is given by 

V{<i),x) =^ 1^ {^oA + [^Jal + \^Hl + XH^ e"^"^ - + ^2^2 ^-Axj 

- (^^Jal + X^Hl + \H^ e-^^ + + A^ff^ g-^xj ^ (45) 
For the model (41) with (34), the eigenvalues (14) are given by 



A H Ae~^* 

M^ = M^ = -3+-, My = -3+^j^. (46) 



If 3 > and 3 > ^^'^ — , the solution is stable. 

Since the solution is unstable if 3 < or 3 < ^^'^2 — 1 there is a possibility that the universe could evolve to the 
de Sitter space-time as in the model (37). For the potential (45), there is an extremum when 



,1 + X^Hl + \H^ e^^ + ^al + X^Hfe^^ = + XHo , (47) 

where the value of the potential V{^, x) is identical with that in (40) and positive if 12iJo > A. The solution (47) is 
not stable, again. 

As one more example, we consider the realistic model which contains the inflation at f ^ — 00, phantom crossing 
at t = 0, and the Little Rip when t — )• 00: 

H = HocoshXt. (48) 

Here Hq and A are positive constants. Since 

H = HoXsinhXt, (49) 
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we find H < when t < 0, that is, the universe is in non-phantom phase and H > when t > 0, that is, the universe 
is in phantom phase. There occurs the phantom crossing at t = 0. Therefore the present universe corresponds to 
t 0. When At 1, we find that the Hubble rate H behaves as 

H ^ f^e^* , (50) 

and therefore there occurs the Little Rip. The EoS parameter w is now given by 

2AsinhAi , , 

w = -1 5 — . (51) 

2Ho cosh^ Xt ^ 

Hence, w < —1 when t > and w > —1 when t < 0. In the limit t — > ±oo, w — > —1. Thus when t — > — oo, there 
occurs the accelerating expansion, which may correspond to the inflation in the early universe. When w = that 
is, 

AsinhAf ^ , „, 

2—=-!' (52) 

Hq cosh Xt 

there occurs the transition between non-accelerating expansion and accelerating expansion. There are two negative 

solutions in (52) in general. Let us denote the solution as ti and t; and assume ti < ti < 0. Then t = ti corresponds to 
the end of inflation and t = ti to the transition from the non-accelerating expansion to the late accelerating expansion 
in the present universe. More explicitly 



A // A \^ . . , , A 




Let the present universe corresponds to t = ipresent- Since ti — tpiosont = 137 x 10^ ~ -jj-^ — = 138 x 10* years 

^ -"present 

(^present ~ 70km/s Mpc) , if we assume tpresent = 0, however, Eq. (52) does not have a solution for fj. Then we may 
assume ^present > 0, that is, the present universe is after the phantom crossing. In this case in principle, one can solve 
(52) with respect to A. Then we can obtain the value of ti from the second equation in (53). Roughly one can expect 
the magnitude of the value could be 50 x 10* years and therefore the realistic cosmology follows. 
The acceleration by the gravitational force between the sun and the earth is given by 

ag = . (54) 

Here I is the distance between the sun and the earth and oja is the angular speed 

27r 



1 year 



= 1.99x10"^-^ (55) 



If the acceleration of the inertial force by the expansion (32) exceeds Og, there occurs the rip between the earth 
and the sun, that is 



IH^ ~ ^e^^* > ag . (56) 



which tells e^-^* = 8.35 x 10^^ or 



Xt = 22.9 . (57) 



If A = O (10-1") years-\ t ~ lO" years. 

For the model (48), by choosing a in (5) to be a constant a = HqX, one finds 



2HoX ^ ^ 2iJoA 

Using (6), it follows 



.;(</>) = ^e-^^ v{x) = ^ cosh Ax. (58) 



f{ct>, x) = Ho {e-^* + sinh Ax) , (59) 
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TABLE I: The (in)stability of the solutions in the models 



Models 


Stability of the reconstructed solution 


Existence of de Sitter solution 


Stability of de Sitter solution 


(21) 


stable 


no 




(35) 


stable if 3 > ^ 
^ xir- 


yes if 12 > ^ 


unstable 


(43) 


stablc if 3 > A and 3 > -^^^ — 


yes if 12 > ^ 


unstable 


(58) 


stable if 3 > ^ 


no 





and the potential in (8) is given by 

^i^'X) = ^ {h^ {e-^^ + sinhAx)' - iToAe"^* + HoXcoshXx] . (60) 
For the model (50) with a = HqX, the eigenvalues (14) are given by 



Therefore if 3 > all the eigenvalues are negative and therefore the solution is stable. 

For the potential (60), there is no extremum and therefore there does not exist the solution corresponding to the 
de Sitter space-time. 

Thus, we constructed scalar models which describe the cosmological solutions with and without Little Rip and 
investigated the (in)stability of the solutions. We also investigated the existence of the solution describing de Sitter 
space-time and the stability of the de Sitter solution when it exists as well as possible transition of Little Rip cosmology 
to de Sitter one. The results are summarized in Table L 



IV. RECONSTRUCTION IN TERMS OF E-FOLDINGS AND SOLUTION FLOW 



Here we consider the reconstruction of the two scalar model in terms of the e- foldings N. (For such a formulation 
in modified gravity, see [15].) We also investigate the flow of the solution by defining dimensionless variables, which 
give the fixed points for some solutions. 

Let us consider two scalar model again. By using the e-foldings N the FRW equations (2) and the scalar field 
equations (9) are rewritten as 



1 



-w(0)(/) +-^ri{x)x + 



= uj{(t)) 

= ??(x) 



y + 3 



ir_ 



x"+(3+f)x']+^.,x(x)x'^ + ^ 
where ' denotes the derivative with respect to the e-foldings = Ina. New function f{(l>,x) is defined by 

d(t)uj{(l)) + j dx??(x) 



f{<t>,x) = /oexp<{ -- 



Here /o is a dimensionless constant. We assume the potential V{(j),x) is given by 



mx) = 



3/(</>,x)' 



(62) 
(63) 
(64) 

(65) 
(66) 
(67) 
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Let 



~ / N N 

= (68) 



Then if the functions w(^) and ry(x) satisiy the following relations 



a solution of 4>, x and H is given by 



K K 



Then one can obtain a model which reproduces arbitrary expansion history of the universe given hy H = f{N)/K, by 
choosing ui{(j)), t]{x) and V{(j),x) by (67) and (69). 

Note that the solution (70) is one of the solutions in the model (1) with (67) and (69). In order to consider the 
structure of the space of the solutions, by defining dimensionless variables, we investigate the flow of the general 
solutions. Besides the solution (70), in general, there are other solutions including the one describing the de Sitter 
space corresponding to the extrema of the potential V{(t), x), where cj) and x are constant. We choose the variables so 
that both of the solution (70) and the de Sitter solution correspond to fixed points. 

Let us introduce the dimensionless variables as follows: 

X = K<^', Y = Kx', ^=-, W = k{cI)-x)- (71) 

X 

Eqs. (64) and (65) are rewritten as 



X' = 3{l-X) + -io{cf>{Z, W)) {X - X^) + -Mx{Z, W)) {X - Y') 
u:,^{4>{Z, W)) 6 (1 - X^) - rjixiZ, W)) (V^ - X') 



^2KUj{<t>{Z, W)) 6 - oj{<j){Z, W)) - vixiZ, W)) 
Y' = 3 {!-¥) + ^ooi^iZ, W)) {¥ - X^) + ^vixiZ, W)) {Y - Y^) 
^ tjMZ, W)) 6 (1 - Y^) - uj{4>{Z, W)) {X^ - Y^ 



(72) 



(73) 



2n TjixiZ, W)) 6 - a;(<A(Z, W)) - r?(x(Z, W)) 

^,_ ix-yzw-z) ^ (74, 

W' = X-Y. (75) 



Now the Hubble rate H is given by 



f{cl,{Z,W),x{Z,W)) / (!>-u:{c^{Z,W))-ii{x{Z,W)) 

K ^ 6-u{(l){Z,W))X^-jj{x{Z,W))Y^' ^ ' 

In order for the Hubble rate to be real, the values of X,Y, Z and W are restricted to a region 

6-u;{^{Z,W))X^-r,{x{Z,W))Y^ 



6-uj{cl>{Z,W))-v{x{Z,W)) 

When a;(^) and r/(x) satisfy (69), this system has two fixed points as follows: 

Point A {X, Y, Z, W) = (1, 1, 1, 0) 

Here the solution is given by (70). 

Point B {X,Y,Z,W) = (/3i,/3i,l,0) 
We now define (3{N) by 



> . (77) 



(3{N) = Po + f3iN . 



(78) 
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Here /3o and Pi are dimensionless constants which satisfy the identities 



6 



mm)). 



3f{l3{N)f\ 



kuj{P{N)/k) k7?(/3(7V)/k) 1 + /3i 
This point exists only if there exist /3o and /3i which satisfy (79). In this point, the solution is given by 



= X = 



H 



f{P{N)) SmN)) + f'iP{N)) 



K Yww)+AV'(/3(iv))' 

Especially when (5\ = 0, this point describes de Sitter space-time. 
We now choose w(^) and ■q{x) as 



(79) 



(80) 



^{(P) T > U , 

/(k0) 

, , f(/^x) + Vf(Kx)^ + a(^x)^ 
= 70^^^^ 



(81) 

(82) 



where a{N) is an arbitrary function. If we choose a{N) = aof'{N) and f{N) being a monotonically increasing or 
decreasing function, co{<j)) and ry(x) are respectively given by 



;(^) = - (1 - e^/l + 



vix) = - (1 + 



f{^x) 



(83) 



Here ao is a dimensionless constant and e = f'{N)/\f\N)\. Then l^(?ii, x) has the following form: 



v{<t>,x) = -f{^^<t>)^-''^^f{i^x) 



l+eVl+ao^ 



1 - eVl + QQ^ f («<^) 1 + eVl + ao' fjKx) 
6 /(k0) + 6 fiKx) 



(84) 



As an example, we consider 

f{N)^foe^^, (85) 
where A is a dimensionless constant. Eq. (84) has a solution H{N) = foe^^ /k. Then w((/)) and ri{x) are given by 

a;(</.) = - (l-eVl + ao') A, rj{x) = - (l + e^^ 1 + ao^) X , (86) 

and y(0, x) is: 



^^(0,X) 

The EoS parameter of the point A is 



^ (^1 + ^) exp [(1 - e^l + c^o') Ak<^ + (l + e^l + "o') Akx 



«^ = -l-3A, 



(87) 



(88) 



which is independent of A''. If A > 0, the point B exists and is located in (3/A, 3/A, 1,0). Then the solution is given 
by 



X = 



3{N_No}_ 
Xk 



H 



/o./A 



,3(iV-JVo) 



K V 3 



(89) 



Here No is an arbitrary constant. The EoS parameter of this point is —3. In this model, the dynamics of X and Y 
are independent of Z and W. Therefore, we consider a small fluctuation from each fixed point by 



X{N)=Xo + dX{N), Y{N)=Yo + 6Y{N). 



(90) 
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FIG. 1: Each vector denotes (X'/50, F'/SO), which is independent of Z and VK. The parameters are A = 3/2 and ao = 1. The 
point A is located in (1, 1), where the EoS parameter is —2. The point B is located in (2,2), where the EoS parameter is —3. 
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FIG. 2: Each vector denotes (X'/20, Y'/20), which is independent of Z and W. The parameters are A = —3/2 and ao = 1. 
The point A is located in (1, 1), where the EoS parameter is 0. 



Here Xq and Yb are the values of X and Y in each fixed point. Then (72) and (73) have the following form: 



8X' 
. SY' 



-3 - A + (1 - eVl + ao^) AXq (l + ey/l + ao'^) AYq 



-3 - A + (1 + eVl + ao') XYq I \6Y 



(91) 



\^ (1 - eVl + ao^) XXo 
The eigenvalues of this matrix (91) are given by 

c7i = -3-A, (72 = -3 - A + (l - eymv) AXo + (l + eyr4W) AFo • (92) 

These indicate that the point A is stable if —3 < A < 3 and unstable if A < —3 or A > 3. Similarly, the point B is 
stable if A > 3 and unstable if < A < 3. The dynamics of X and Y are shown in FIG. 1 and FIG. 2. 
We now consider another example 



f{N) = foN^.^ 



(93) 



where 7 is a dimensionless constant. We should note that the case 7=1 corresponds to the Little Rip model (17). 
Then Lo{(j)) and r]{x) are given by 



a;(<^) = - (1 - €^/l + ao^) ^ , vix) = - (l + eVl + ao^) ^ , 
and V{(j),x) has the following form: 



(94) 



3/0^ 



l-eVl+ao^ 



(«x) 



1 + 



1 — e\/l + ao^ 7 



+ 



1 + eVl + ao^ 7 
6 Kx. 



(95) 
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FIG. 3: Each vector denotes (X'/50, K'/SO) with Z = \ and VK = = x = 1/^)- The parameters are 7 = 1 and qo = 1. 
The point A is located in (1, 1), which corresponds to the Little Rip universe. 



The EoS parameter of the point A is 



- = -l-|^, (96) 



which becomes —1 when N ^ 00. If 7 < 1/2, the point B exists and is located in (0, 0, 1, 0). Then the solution is 

, = ,= 1^, H^J^(^]"'-\ (97) 



6k ' V6k\1-2^ 
This point corresponds to the de Sitter space-time. We consider a small fluctuation from the point B by 

</'(^)-m^ + '5</.(iV), x{N) = ^^^ + dxiN). (98) 



Then (72) and (73) have the following form: 

= fl-27)M +JsA_is(^-(^-''{II^'^ -{1 + eVT+a^A A^' 

V<^^7 W) V - (1 - 7 1 - (1 + eVTT^ 7 J [sXj 

The solution of this equation (99) is given by 



(99) 



(^A = ( 1 + eV l + gp^ ^ Ci+e'-^+^ + Ci-e"^-^ /l\ C^+e'^+^ + Cz- e"^-^ 

\Sx) \-l + eVT+^) « VV K ■ ^ ' 

Here Ci± and C2± are arbitrary constants and ai± and cr2± are given by 

-3±3v/l + 8(l-27) -3±3Vl + 8(l-27)2 

^^± = — — ' """^^ — — ■ ^ ^ 

These indicate that the point B is always unstable. On the other hand, when ^, x — > 00, the solution of X and Y are 
given by 

a: = 1 - ATi , Y = l^Yi e^^^ . (102) 

Here Xi and Yi are arbitrary constants. This indicates that the values of X and Y approach to (1, 1) if ^, x — 00 
when N ^ 00. Then the Hubble rate approaches to /qN'^/k, which corresponds to the Little Rip universe if 7 = 1. 
The dynamics of X and Y arc shown in FIG. 3 and FIG. 4. 

Then we have completed the formulation of the reconstruction in terms of the e-foldings. The e-foldings description 
is directly related with redshift and therefore, with the cosmological observations. We also investigated the flow for the 
solution, which shows the (in)stability of the reconstructed solution obtained in a large range. Even if the solution is 
stable, when the stable region is small, the evolution of the universe depends strongly on the initial conditions. If the 
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FIG. 4: Each vector denotes {X' /2Q, Y' /2Q) with Z = 1 and W = (c^ = X = 1/2k). The parameters are 7 = -1 and ao = 1. 
The point A is located in (1, 1). The point B is located in (0,0), which corresponds to the de Sitter universe. 
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FIG. 5: Each vector denotes (Z'/20, W' which is independent of the form of w(0), ri{x) and V(0, x)- The dynamics of Z 
and W are classified into four types according to the values of X and y as (a) < y, X < y, (b) < y < X, (c) > y, 
X > y and (d) > y > X. The fixed points are located in (1, 0). 



initial condition is out of the range, the imiverse does not always evolve to the sohition obtained by the reconstruction. 
On the other hand, if the stable region is large enough, even if the universe started from an initial condition in a 
rather large region, the universe evolves to the solution obtained by the reconstruction. 

As mentioned above, the model (93) with 7 = 1 corresponds to the Little Rip model (17). In fact, Eq. (17) indicates 
that 



(103) 
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and therefore 

H = XN. (104) 

By comparing (104) with (93), we find /o = A. Then as written after (19), the parameter /o is bounded as 2.74 x 

10~^ Gyr~^ < 2/o/\/3 < 9.67 x 10~^ Gyr~^ by the results of the Supernova Cosmology Project [16]. If we choose so 
that the present universe corresponds to N = Nq, we have foNo ~ 70km/sMpc. 

V. DISCUSSION 

In summary, we gave a general formulation of reconstruction in two scalar model and investigated the stability of 
the solution. This formulation helps us to construct a model which has a stable cosmological solution describing the 
phantom-divide crossing. By using the formulation, we constructed a model which describe the cosmological solutions 
with and without Little Rip and investigated the (in)stability of the solutions. The existence of the solution describing 
de Sitter space-time was also investigated and furthermore the stability of the de Sitter solution when it exists as well 
as possible transition of Little Rip cosmology to dc Sitter one was investigated. We also considered the reconstruction 
of the two scalar model in terms of the e- foldings N and investigated the flow of the solution by defining dimensionless 
variables, which give the fixed points for some solutions. 

Finally, let us make several remarks about the relation of the qualitative behavior of the Universe evolution and 
the shape of the scalar potential. 

In case of the usual canonical scalar field as <j> in (1), when the field climbs up the potential, the kinetic energy 
decreases until the kinetic energy vanishes. Even in case of the phantom field as x in (1) with non-canonical kinetic 
term, the kinetic energy decreases when the field climbs up the potential. In case of the phantom field, the kinetic 
energy is unbounded below and therefore the absolute value of the kinetic energy increases when the field climbs up 
the potential. The big rip or Little Rip occurs when the potential goes to infinity. If the potential tends to infinity 
in the finite future, the evolution corresponds to the big rip but if the potential goes to infinity in the infinite future, 
the evolution corresponds to the Little Rip. Then the necessary condition that the big or Little Rip could occur is 

1. The potential does not have maximum and it goes to infinity. 

2. There is a path in the potential that the potential becomes infinite but the kinetic energy of the canonical scalar 
field is vanishing. 

Since we identify the scalar field (j) and x with the cosmological time, the second condition means w(^) in (1) goes 
to zero when (p goes to infinity and therefore the phantom field x dominates. Conversely, if there is a maximum in 
the potential or there is no path in the potential that the potential becomes infinite but the kinetic energy of the 
canonical scalar field goes to zero, there does not occur big rip nor Little Rip. 

Let us suppose the case that there is a maximum in the potential. If the fields stay near the potential maximum, 
the universe becomes asymptotically de Sitter space-time. If the fields go through the maximum and the potential 
decreases, the kinetic energy of the canonical scalar field increases but the absolute value of the kinetic energy of the 
phantom field decreases. If the kinetic energy of the phantom field goes to zero, the canonical field becomes dominant 
and the Universe could enter the non-phantom (quintessence) phase and there might occur the deceleration phase in 
future. 



Acknowledgments 

We are grateful to M. Chaichian, P. Frampton, V. Mukhanov, M. Sasaki, and R. Scherrer for related discussions. 
This research has been supported in part by MEC (Spain) project FIS2006-02842 and AGAUR(Catalonia) 2009SGR- 

994 (SDO), by Global COE Program of Nagoya University (G07) provided by the Ministry of Education, Culture, 
Sports, Science & Technology and by the JSPS Grant-in-Aid for Scientific Research (S) # 22224003 and (C) # 
23540296 (SN). 



[1] R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 071301 (2003) [astro-ph/0302506] . 
[2] S. 'i. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71, 063004 (2005) [hcp-th/0501025]. 
[3] S. 'i. Nojiri and S. D. Odintsov, Phys. Rapt. 505, 59 (2011) [arXiv:1011.0544 [gr-qc]]. 



15 



[4] P. H. Frampton, K. J. Ludwick, R. J. Schcrrcr, [arXiv:1106.4996 [astro-ph.CO]]. 

[5] P. H. Frampton, K. J. Ludwick, S. Nojiri, S. D. Odintsov and R. J. Schcrrcr, arXiv:1108.0067 [hcp-th]. 
[6] S. Nojiri, S. D. Odintsov and D. Saoz-Gomcz, arXiv:1108.0767 [hcp-th]. 

[7] I. Brevik, E. EUzalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 84, 103508 (2011) [arXiv: 1107.4642 [hep-th]]. 

[8] P. H. Frampton and K. J. Ludwick, Eur. Phys. J. C 71, 1735 (2011) [arXiv:1103.2480 [hep-th]]. 

[9] L. N. Granda and E. Loaiza, arXiv:1111.2454 [hep-th]. 
[10] B. Feng, X. -L. Wang and X. -M. Zhang, Phys. Lett. B 607, 35 (2005) [arXiv:astro-ph/0404224 [astro-ph]]. 
[11] E. Elizalde, S. 'i. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 043539 (2004) [hep-th/0405034]. 
[12] Y. -F. Cai, E. N. Saridakis, M. R. Sctaro and J. -Q. Xia, Phys. Rcpt. 493, 1 (2010) [arXiv:0909.2776 [hep-th]]. 
[13] S. Nojiri and S. D. Odintsov, Gen. Rol. Grav. 38, 1285 (2006) [arXiv:hep-th/0506212]. 
[14] S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 632, 597 (2006) [arXiv:hep-th/0507182]. 
[15] S. 'i. Nojiri, S. D. Odintsov and D. Saez-Gomez, Phys. Lett. B 681, 74 (2009) [arXiv:0908.1269 [hep-th]]. 
[16] R. Amanullah, C. Lidman, D. Rubin, G. Aldering, P. Astier, K. Barbary, M. S. Burns, A. Conley et al., Astrophys. J. 
716, 712-738 (2010). [arXiv: 1004. 1711 [astro-ph.CO]]. 



